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Abstract. The physics of energy transfer between the laser and the plasma in laser 
wakefield accelerators is studied. We find that wake excitation by arbitrary laser shapes 
can be parameterized using the total pulse energy and pulse depletion length. A tech- 
nique for determining laser profiles that produce the required plasma excitation is 
developed. We show that by properly shaping the longitudinal profile of the driving 
laser pulse, it is possible to maximize both the transformer ratio and the wake ampli- 
tude, achieving optimal laser-plasma coupling. The corresponding family of laser pulse 
shapes is derived in the nonlinear regime of laser-plasma interaction. Such shapes 
provide theoretical upper limit on the magnitude of the wakefield and efficiency of 
the accelerating stage by allowing for uniform photon deceleration inside the laser 
pulse. We also construct realistic optimal pulse shapes that can be produced in finite- 
bandwidth laser systems and propose a two-pulse wake amplification scheme using the 
optimal solution. 



INTRODUCTION 

Recent advances in laser technology allow one to create laser pulses with virtually 
arbitrary temporal intensity profiles using amplitude and phase shapers [1-3]. Such 
laser pulses with non-Gaussian axial intensity are now being considered for applica- 
tions as drivers in Laser Wakefield Accelerators (LWFA). Shaped lasers provide the 
means of controlling the generation of plasma wake and thus offer the possibility 
of optimization of wake excitation and accelerating efficiency However, progress 
in finding "the optimal" shape has been hindered by the apparent complexity of 
the problem. Not only is the parameter space of possible shape functions huge, 
but also the generated wakefield is a nonlinear function of laser intensity, requir- 
ing numerical solution of differential equations in a variational calculation. As a 
result, several groups turned to trial and error methods such as genetic algorithms 
for optimization [2,3]. Still, even these methods require consistent classification 
of laser shapes so that different pulses can be meaningfully cross-compared while 
desired properties such as wake amplitude or efficiency are optimized. In this paper 
we reanalyze the process of wake generation and argue that the only two physical 



parameters that describe a laser shape from the stand point of wake excitation are 
the total pulse energy and its depletion length. Using these parameters we find the 
analytic expression for the family of optimal laser shapes that maximize both the 
wakefield and the accelerating efficiency We also develop a method for determining 
the shape of a laser that produces a required value of wakefield without explicitly 
solving the wake equation. This opens the way for obtaining laser shapes that 
satisfy other optimization criteria specific to given experimental conditions. 



Wakefield accelerators such as the laser-driven LWFA [4] or electron beam driven 
Plasma Wakefield Accelerator (PWFA) [5] can be viewed as two-step energy trans- 
fer systems: in the first step the driver deposits energy into wake excitation of the 
plasma, and in the second step the energy is taken from the wake by the accel- 
erating beam. While the second step is the same for both accelerating schemes, 
the physics of driver energy deposition is quite different between them. In PWFA 
the electron beam loses energy to the plasma through interaction with the in- 
duced electrostatic field, while in the LWFA laser energy loss occurs via photon 
red-shift or deceleration [6]. This process can be understood as follows. Podero- 
motive force of the laser modifies both the density n e and the Lorentz factor 7 
of the plasma electrons. This produces modulations in the nonlinear index of the 



refraction 77 = [1 — (u p / 'u) 2 n e j '^n p } 1 / 2 , where cu p = ^Aive 2 n p /m e is the unper- 

turbed plasma frequency, and uj is the frequency of the laser. The wake-induced 
modulations of the refractive index appear stationary in the reference frame co- 
moving with the laser and cause laser photons to red- or blue-shift depending on 
the sign of refractive index gradient [7]. Due to negligible scattering in the set- 
ting of laser accelerators, the photon number in the laser is essentially constant, 
and the energy deposition into the plasma is therefore determined by the pho- 
ton deceleration. To address this quantitatively we consider a laser propagating 
along the z axis with initial frequency u ^> uj p . In the laser comoving frame, the 
plasma response can be written in terms of the independent dimensionless variables 
( = k p (z — v g i) and r = k p ct, where k p is the plasma wavenumber, and v g ~ — c 
is the laser group velocity (for convenience, the laser is moving in the negative z 
direction). Introducing dimensionless normalized scalar and vector potentials <p(() 
and a((), the parallel and perpendicular electric fields are E\\ = — (mc 2 k p /e)d(f)/d( 
and E± = —(mc/e)da/dt = —(mc 2 k p /e)da/d(. The wakefield generation equation 
can then be written as [8,9]: 



ENERGY TRANSFER IN LWFA 





(1) 



where x = 1+0 is the modified electrostatic potential, and a 2 (C) is the dimensionless 
laser intensity envelope averaged over fast oscillations. Prior to the arrival of the 



laser the normalized wakefield £ = eE^/maj p = —dx/d( is zero. A formal solution 
for the electric field outside the laser can be written as the first integral of (1): 
[£ out (()] 2 = — (x — l) 2 /x+ Jf^ a 2 x' /x 2 d(, which reaches a maximum value at x — 1: 

K^ = -£V(c)(|iK (2) 

This expression can be understood in terms of the deposition of laser energy into 
plasma. For this we use the formula for local frequency shift of laser photons 
obtained from the analysis of laser evolution equation [10,11]: 

The energy density in the wake from (2) can then be interpreted as the intensity- 
weighted integral of the photon deceleration throughout the pulse. Let's denote 
the wake-dependent part of the photon deceleration function as k(() = x' /x 2 . The 
value of the peak wakefield in (2) is then bounded from above by the total laser 
energy (the integral of a 2 ) and the maximum photon deceleration n max : 

/OO /'OO 
a 2 {C)K{QdC<K max a 2 (C)dC (4) 
-oo J —oo 

where n max is the maximum of k(Q inside the laser. Maximum photon decelera- 
tion n max actually has a simple physical interpretation. It is closely related to the 
characteristic laser depletion length l d , or the distance in which the maximally de- 
celerated laser slice red-shifts down to uj v (assuming no evolution of the wakefield). 
From (3) this characteristic depletion length is: 

Id = [(uJo/ujp) 2 - l]/k p K max . (5) 

The peak wakefield outside the laser then scales with depletion length l d and di- 
mensionless pulse energy e = Jf^ a 2 (()d( as: £°^ x < ^[{uiq/uj p ) 2 — l]e /k p l d . The 
range of achievable wakefields is therefore set by the total pulse energy and its 
depletion length. For pulses of fixed energy and depletion length the actual value 
of the wakefield within this range will depend only on particular laser shape, and 
can be optimized by varying the shape subject to constraints. 



WAKEFIELD OPTIMIZATION 

One possible optimization problem can be formulated as follows: given some 
fixed laser energy e and the desired depletion length (K max = k ), what laser 
shape would produce the largest possible wakefield? According to (4) this amounts 
to finding a shape that maintains the largest for the duration of the pulse. 
Since the maximum k(Q is fixed by the depletion length, such a shape should have 
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FIGURE 1. General shape of the nonlinear optimal laser intensity profile and its corresponding 
wakefield (arbitrary units) 

a constant photon deceleration throughout the pulse, k(() = k . If the laser is 
present for ( > 0, then in order to satisfy the boundary condition of quiescent 
plasma before the laser, the photon deceleration should rise from to value at 
the very beginning of the pulse, e.g., like a step-function: = /to#(C + )- Here, 
( + = ( — + in order to avoid ambiguities with the values of step-function at 0. 
The corresponding laser profile is then found from the wake equation (1): 

_ 2^( C +) 4«gfl2(C+) 1 
a/(C) "(l--oC) 4 + (l--oC) 5+ (l--oC) 2 ' (6) 

where £ £ [0, C/ < V^o], an d ^(C + ) is a delta-function such that Jq >0 S(y + )dy = 1. 
A schematic drawing of the optimal laser intensity variation and its associated 
plasma wakefield are shown in Fig. 1. Generally, the shape consists of a 5-function 
at the front followed by a ramp in intensity which is cut off at Q. In the linear 
regime, when a 2 « 1, Ko -* 0, the ramp reduces to a triangular shape found 
in [11,12]: a 2 = 2k (5(( + ) + (). We note that (6) describes a family of shapes, 
rather than a fixed shape. The actual profile of the optimal pulse depends on the 
deceleration parameter kq set by the desired depletion length and the pulse length 
(f, which is determined from the available total energy: 

Although the pulse length cannot exceed ( c = 1/k,q, the rise of a 2 towards the 
end of the pulse guarantees that any finite laser energy can be accommodated for 
Q < ( c . The two terms in (7) represent the energy contained in the 5-function 
precursor and the main pulse. It is clear that for a fixed total energy there exists 
a maximum value of k = e /2 which is achieved when (f — > 0, i.e., all of the 
energy is concentrated in the 5-function. This shape, which is a particular case of 
the general optimal shape (6), excites the largest possible wakefield and has the 



smallest depletion length among all pulses of fixed energy. For circularly polarized 
pulses with cylindrical transverse crossection of radius r and wavelength A, the 
maximum achievable wake is then given by: 



E max = 6.5AE wb 
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where Uq is the total pulse energy (in Joules) and E w b = 96[n p /10 18 cm 3 ]GV/m is 
the nonrelativistic wavebreaking field. 



EFFICIENCY OPTIMIZATION 

While generation of large accelerating gradients is a prerequisite for a successful 
accelerating scheme, the efficiency of acceleration should also be considered. For 
an accelerating scheme that involves transfer of energy from the driver beam to 
the accelerating beam, efficiency is measured in terms of the transformer ratio, or 
the ratio of the maximum rate of energy gain per particle of accelerating beam 
to the maximum rate of energy loss per particle of the driving beam. In the case 
of laser-plasma accelerators, where the driving and accelerating beams consist of 
particles of different species, the following kinematic definition is more useful: 

o _ \dla/dz\ 

max /„\ 

\oid/oz\ max 

where 7^ and 7 a are Lorentz factors for the driving and accelerating beams. In 
LWFA the particles in the trailing electron bunch are accelerated via electrostatic 
interaction with the wake, so \dj a / dz\ max = \eE™ ax \lm e (? = k p \£™* x \. For the 
laser propagating in plasma 7^ uj/u p , so \d^d/dz\ is the photon frequency shift 
given by (3). The transformer ratio for LWFA is then: 

dLWFA _ 2u; \d x /dC\max ^ \ c\out , , / ln \ 
U - l^\in ^ \ C \max K P l d- l iU J 



cu p \d(l/x)/d(\ 



max 



Defined this way the transformer ratio can have several interpretations. On the 
one hand, it is a measure of local accelerating efficiency, or the amount of increase 
in 7 of the accelerating electron per unit loss of 7 of the laser. On the other hand, 
transformer ratio is proportional to the maximum energy that can be transferred 
to the accelerating beam particle over the laser depletion length (assuming no 
evolution of the wake during propagation). Therefore, an efficient accelerating 
scheme should attempt to maximize the transformer ratio. 

There are several ways to find the laser shape that maximizes R. Among the 
pulses of fixed energy and depletion length, R is maximized by a pulse that produces 
the largest wakefield as can be seen from (10). The optimal shape found in Eq. (6) 
satisfies this requirement. Alternatively, one can relax the energy constraint and 



instead look for a laser profile that has the largest depletion length among all the 
shapes that produce a given maximum wakefield behind the laser. Although this 
reasoning leads to the same resulting shape, we include the proof for completeness 
as it demonstrates a useful technique for determining laser shapes that satisfy 
constraints on the values of the wakefield. 

In order to find the shape that maximizes the transformer ratio, we vary the 
photon deceleration function k,(() inside the laser. We require that k,(() be positive 
definite, i.e., laser photons only lose energy to the plasma and do not reabsorb 
energy from the wake. The advantage of varying k(Q rather than a 2 (() directly 
is that one can immediately write down the solution for the wakefield potential 
x(() in terms of the photon energy deposition function = /i^ K,((i)d(i, i.e., 
x(Q = 1/(1 — fp(C))- The corresponding laser shape is then determined from the 
wakefield equation (1): 

„.(0 = (2,»(0 + 1*0' - i = ^ + ^% + ^ - MID 

Note that not all functions V'(C) should produce physical, i.e., positive a 2 ((), 
and the validity of a given k(Q should be checked through (11). By considering 
photon energy deposition in the pulse all possible laser shapes that produce a 
given wakefield can be mapped onto a bounded space of monotonically increasing 
functions ip(C), whose end values on the interval [0, Q] and derivatives at Q are 
constrained by the required maximum value of the wakefield. From the first integral 
of the wakefield equation we can relate the modified potential Xf = x(Q) and the 
electric field x'f = x'(Q) evaluated at the end of the pulse: 

Monotonicity of V'(C) follows from the requirement k{Q > 0, and the bounds 
on V'(C) are = V'(O) — V'(C) — ^max < 1. A few sample solutions for V'(C) an d 
corresponding photon deceleration and laser shapes are plotted in Fig. 2. The func- 
tion that results in the largest transformer ratio should possess the smallest 
maximal slope in the interval [0, Q] - this will maximize the depletion length for 
a fixed £™l x - Such curve is unique and is represented by curve 2 in figure 2. It 
is a straight line with slope ip'(Cf) = VKC/Ocrit/C/ > where the value of V'(C/)crit is 
determined from substituting ip'((f) into eq. (12). Let's show that this line has the 
smallest maximum slope. Since ip'(Cf) is a decreasing function of ip{Cf) ior a fixed 
£max ( ec l- (12)), all curves V'(C) with ip(Cf) < V , (C/)ca-it (such as curve 1 in figure 2) 
will automatically have larger slope at Q: > V'(C/)crit/C/- O n the other hand, 

the curves with ip((f) > V'(C/)crit (such as curves 3 and 4) should have a slope larger 
than ■0(C/)crit somewhere between and Q in order to be larger than ■0(C/)crit & 1 
(f. We therefore prove that the function = Kq(, where Kq = V 7 (C/)crit/C/i is an 
integrated photon deceleration profile that maximizes the transformer ratio. 
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FIGURE 2. a) Sample photon energy deposition ^>(£) for pulses of length not exceeding Q = 2n: 
b) corresponding photon deceleration functions; c) resulting laser intensity profiles. All shapes 
produce the same maximum wakcficld. 



The photon deceleration function associated with this V'(C) i s a constant /c(£) = 
kq and the resulting laser shape is the same as given by (6). The optimal trans- 
former ratio associated with this shape can be found from (10): 



i+(k p L P ni-K (k p L p ]i^ 



jj>LWFA _ ^ 

uj p \^ [1 - K (kpL p )] 



where L p = Q/k p is the pulse length. In the linear regime optimal transformer 
ratios for both LWFA and PWFA schemes scale identically with the pulse/beam 

length: R LWFA - (2u/u p )yfl + {k p L pulse )\ i? PWFA - + {k p L beam f [13]. The 
LWFA scheme is intrinsically more efficient by a factor of 2u/u p , which is needed for 
viability of LWFA since lasers are typically "slower" drivers than electron beams. 



UTILITY OF PULSE SHAPING 

The advantage of using the optimal pulse shape is best seen in comparison with 
the unshaped (Gaussian) pulse. For a given Gaussian pulse (or any other non- 
optimal shape) one can always construct a corresponding optimally shaped pulse 
with the same laser energy such that the photon deceleration across the optimal 
pulse equals to the peak photon deceleration in the unshaped one (i.e., both pulses 
have equal depletion lengths). Unshaped pulses deplete first in the region where 
photon deceleration is the largest, whereas a laser with the optimal shape loses all 
its energy in a depletion length due to uniform photon deceleration, thus enhanc- 
ing instantaneous energy deposition and wakefield. For a numerical example, we 
consider the optimal and Gaussian pulses of total energy 0.5J, wavelength 1/im and 
cylindrical radius 10/im in a plasma with n p = 10 18 cm~ 3 . The transformer ratio, 
the maximum wakefield, the required pulse length, and the corresponding peak do 
are shown in Fig. 3 as a function of depletion length. 
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FIGURE 3. Comparison of the transformer ratio, maximum wakefield, pulse length, and maxi- 
mum normalized vector potential in shaped (diamonds) and Gaussian (triangles) pulses of equal 
depletion lengths and constant pulse energy of 0.5J. 



From Fig. 3 we see that the transformer ratio and the maximum wakefield are 
consistently larger for shaped pulses. In fact, the lines for optimal pulse wakefield 
and transformer ratio represent the theoretical upper limits for all pulses of given 
energy. The Gaussian pulse achieves a maximum transformer ratio when its length 
(measured here as FWHM) equals 1/2 of the relativistic plasma wavelength. The 
effects of shaping are especially prominent for longer pulses, where Gaussian pulse 
yields almost no wake excitation due to plasma oscillations inside the pulse that 
cause part of the laser photons to absorb energy from the wake. On the other 
hand, a shaped laser postpones plasma oscillation until the end of the pulse, and 
all photons decelerate uniformly. For very short pulses, the differences between the 
two shapes are minimal. This is due to the fact that very short Gaussian pulses of 
fixed energy asymptotically approach the delta-function limit of the short optimal 
shape. For these short pulses the wakefield attains the maximum value given by 
(8) as the depletion length reaches the minimal value for given pulse energy. 

Although short pulses generally produce the largest wakefields, their efficiency is 
close to minimal possible, as the depletion length decreases faster than increase in 
the wake. Therefore, the choice of the appropriate pulse shape for LWFA stage will 
depend on specific experimental conditions. If the laser-plasma interaction distance 
is limited by instabilities, diffraction or dephasing, then in order to maximize the 
electron energy gain one should try to achieve the largest accelerating gradient, 
which can be accomplished with ultrashort pulses. For some regimes of plasma 



density and laser energy available laser systems may be unable to produce pulses 
short enough so that the pump depletion length is longer than the characteristic 
instability or dephasing length. In this case shaping the laser will increase the 
wakefield over the interaction distance, even though it will be below the maximum 
possible if a shorter pulse were used. If the interaction length is less constrained, 
such as the case for propagation in plasma channels [14], then using a finite-length 
shaped pulse will result in a greatly improved overall energy gain per stage as 
can be seen from Fig. 3. An added benefit of pulse shaping is the suppression of 
modulational instability that affects unshaped pulses that are longer than plasma 
wavelength. When all photons red-shift, or "slow down", at the same rate, different 
slices of the laser do not overrun each other, and the ID laser self-modulation is 
suppressed. 



REALISTIC PULSE SHAPING 

As the optimal pulse shape is associated with a delta-function precursor, the 
feasibility of such a structure may be a concern. We note that the purpose of this 
precursor is to bring the photon deceleration from zero in the quiescent plasma 
before the laser to a finite value k,q at the beginning of the main pulse. This can 
also be achieved with a more physical prepulse, whose shape can be found from the 
wake equation once a smooth function k(Q is chosen. 

For our example we choose a photon deceleration function that varies as a 
hyperbolic tangent: k(() = k [1 + tanh(a:(C — Co))]/2, where a is a steepness 
parameter and Co is an arbitrary offset. The photon energy deposition is then 
= Hq[( + ln(cosh(a:(C — ( )))/a}/2, and the corresponding laser shape is found 
from equation (11): 

fl2 /^ = /t asech 2 (a(C - Co)) wg[l + tanh(a(C - Co))] 2 _J 1 ^ 

X 4 (C) ~ " " ' X 5 (C) " X 2 (C) 

where C < C/ an d the function in the denominator is x(C) = 1 + («o/2a) In | + 
Co«o/2 — As before, the pulse length C/ can be found from the total available 

energy of the pulse. By varying a we can change the slope of k(Q as it rises 
from to /to and construct a pulse shape that satisfies experimental constraints 
yet retains essential physics of the optimal shape. For a step-function photon 
deceleration (a — > oo) expression (14) asymptotes to equation (6). However, for 
finite values of a the delta-function precursor spreads out and can even disappear 
as shown in Fig. 4. The family of shapes given by (14) is better suited for the finite- 
bandwidth laser systems that have a lower limit on achievable feature size. The 
values of maximum wakefield for pulses in Fig. 4 is within few percent of the value 
for a delta-function optimal pulse of the same energy and depletion length. This is 
due to the fact that the bulk of the laser pulse still experiences constant maximal 
photon deceleration. The wakefield further degrades with longer rise times of k(C)- 
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FIGURE 4. Laser intensity (shaded) and associated photon deceleration (— k(Q) for pulses of 
the same total energy and characteristic depletion length in the order of increasing a. 

The pulse shaping techniques described so far have assumed that the laser is 
incident on an unperturbed plasma. However, this does not have to be the case, and 
we can construct an optimally-shaped laser that enters the plasma at some phase 
of a pre-existing plasma oscillation. Such oscillation could be left from a precursor 
laser pulse or electron beam as shown in Fig. 5. When there is an existing plasma 
wave, the value x of the modified electrostatic potential at the beginning of the 
optimal pulse will generally be different from unity. In this case the expression for 
the optimal pulse without the delta-function precursor is modified into: 



where we assume that the main pulse lies between (o and Q, and kq = x'((q)/xq. 
If this shape is placed in a correct phase of the oscillation (so that af((o) from 
(15) is positive), it acts as an amplifier of the existing wakefield. The ratio of 
maximum wakefield behind the optimal pulse to the field in front of it scales as 
(R/xq)(uj p /2uj) which for pulse lengths around X p from Fig. 3 can be of order 10. A 
detailed discussion of this scheme and a comparison to the resonant laser-plasma 
accelerator concept [15] will be reported elsewhere. 



As we have shown, the huge phase space involved in shaping laser drivers for 
applications in laser wakefield accelerators can be described using only two param- 
eters: total pulse energy and characteristic depletion length. The shape of photon 
energy deposition (photon deceleration) inside the pulse plays a crucial role for 
both the wake excitation and the evolution of the laser driver. By varying the 
shape of the photon deceleration function for pulses of fixed energy and depletion 
length we were able to optimize both the generated wakefield and the efficiency of 
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FIGURE 5. Laser intensity profiles (a 2 (£), shaded) and normalized electric field for optimally 
shaped main pulse following a Gaussian precursor. 

the accelerating scheme. The method used for obtaining the optimal shapes (6) 
and (14) is actually more general and can be used to determine laser shapes that 
generate other variations in the nonlinear index of refraction. Having a physical 
requirement for the refractive index, which in this case is the requirement of uni- 
formity of photon deceleration, provides a constraint on the functional form of the 
wakefield, which can then be used to find the required laser shape. Alas, such a 
"reverse" solution is not always guaranteed to yield a physical (i.e., positive) a 2 (Q, 
so, in general, caution is advised. 

Several issues should be addressed before the laser pulse shaping concept can 
be fully utilized. Even without the delta-function precursor, the finite laser band- 
width will necessarily smooth out steep rises and falls of the optimal pulse shape. 
Although we do not anticipate adverse effects when the feature size is much smaller 
than the plasma wavelength, the ID self-consistent laser evolution and stability of 
realistic optimal shapes are currently under investigation. Another consideration 
is the influence of the laser-plasma interaction in the transverse dimension on the 
evolution of the pulse. Many of the laser-plasma instabilities are seeded by the 
wakefield-induced perturbations of the index of refraction. As we have demon- 
strated in this paper, the nonlinear index of refraction can be effectively controlled 
through laser shaping, thus suggesting the method of delaying the onset of these in- 
stabilities. Whether this approach increases the growth rates of other instabilities, 
particularly in the transverse dimension, remains to be investigated. 

We would like to thank J. Arons, A. Charman, T. Katsouleas, W. B. Mori, and 
J. Wurtele for fruitful discussions and suggestions. 

REFERENCES 

1. F. Verluise, V. Laude, Z. Cheng, et al, Optics Lett. 25, 575 (2000). 

2. M. Murnane, these proceedings. 

3. M. Downer, these proceedings. 



4. T. Tajima, J. M. Dawson, Phys. Rev. Lett. 43, 267 (1979). 

5. P. Chen, J. M. Dawson, R. Huff, T. Katsouleas, Phys. Rev. Lett. 54, 693 (1985). 

6. S. Wilks, J. M. Dawson, W. B. Mori, T. Katsouleas, M. Jones, Phys. Rev. Lett. 62, 
2600 (1989). 

7. W. B. Mori, IEEE J. Quant. Elec. 33, 1942 (1997) 

8. P. Sprangle, E. Esarey, J. Krall, and G. Joyce, Phys. Rev. Lett. 69, 2200 (1992). 

9. E. Esarey, P. Sprangle, J. Krall, A. Ting, IEEE Trans. Plasma Sci. 24, 252 (1996). 

10. E. Esarey, A. Ting, and P. Sprangle, Phys. Rev. A 42, 3526, (1990). 

11. P. Chen, A. Spitkovsky, AIP Conf. Proc. 472, 321 (1999). 

12. P. Chen, A. Spitkovsky, T. Katsouleas, W. B. Mori, Nucl. Instr. Meth. 410, 488 
(1998). 

13. P. Chen, J. J. Su, J. M. Dawson, K. Bane, and P. Wilson, Phys. Rev. Lett. 56, 1252 
(1986). 

14. E. Esarey, P. Sprangle, J. Krall, A. Ting, G. Joyce, Phys. Fluids B 5, 2690 (1993). 

15. D. Umstadter, E. Esarey, J. Kim, Phys. Rev. Lett 72, 1224 (1994) 



